Analytic bounds on transmission 
probabilities 

Petarpa Boonserm and Matt Visser 

School of Mathematics, Statistics, and Operations Research 
Victoria University of Welhngton, New Zealand 

{petarpa. boonserm, matt. visserjOmsor.vuw.ac.nz 

8 January 2009; ErgX-ed January 7, 2009 

Abstract 

We develop some new analytic bounds on transmission probabilities (and 
the related reflection probabilities and Bogoliubov coefficients) for generic one- 
dimensional scattering problems. To do so we rewrite the Schrodinger equation 
for some complicated potential whose properties we are trying to investigate 
in terms of some simpler potential whose properties are assumed known, plus 
a (possibly large) "shift" in the potential. Doing so permits us to extract con- 
siderable useful information without having to exactly solve the full scattering 
problem. 
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1 Introduction 



In several earlier papers [H [21 El H], the present authors have derived a number 
of rigourous bounds on transmission probabilities (and reflection probabilities, and 
Bogoliubov coefficients) for one-dimensional scattering problems. The derivation of 
these bounds generally proceeds by rewriting the Schrodinger equation in terms of 
some equivalent system of first-order equations, and then analytically bounding the 
growth of certain quantities related to the net flux of particles as one sweeps across 
the potential. In the present article we shall obtain significantly different results, of 
both theoretical and practical interest. 

While a vast amount of effort has gone into studying the Schrodinger equation 



and its scattering properties [5l[6l[3[Hl[9l[ini[Ill[T2l[T3l[Ta[l5l[T6l[I71[ig 



it appears that relatively little work has gone into providing general analytic bounds 
on the transmission probabilities, (as opposed to approximate estimates). The only 
known results as far as we have been able to determine are presented in [11 [21 [3l H] . 
Several quite remarkable bounds were first derived in [1], with further discussion 
and an alternate proof being provided in [2]. These bounds were originally used by 
one of the present authors as a technical step when studying a specific model for 
sonoluminescence [22] , and since then have also been used to place limits on particle 
production in analogue spacetimes [23] and resonant cavities [21], to investigate 
qubit master equations [25], and to motivate further general investigations of one- 
dimensional scattering theory [261 [23, [2H]- Recently, these bounds have also been 
applied to the greybody factors of a Schwarzschild black hole [3]. Most recently, 
significant extensions of the original bounds have been developed [4] by adapting the 
Miller-Good transformations |29j . 

In the current article we again return to this problem, developing a new set of 
techniques that are more amenable to the development of both upper and lower 
bounds on the transmission probabilities. For technical reasons the new techniques 
are also more amenable to investigating behavior "under the barrier". The basic 
idea is to re-cast the Schrodinger equation for some complicated potential whose 
properties we are trying to investigate in terms of some simpler potential whose 
properties are assumed known, plus a (possibly large) "shift" in the potential. 

2 From Schrodinger equation to system of ODEs 

We are interested in the scattering properties of the Schrodinger equation. 




(2.1) 
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where k{x)'^ = 2m[E — V{x)]/h?. As long as V{x) tends to finite (possibly distinct) 
constants V±oo on left and right infinity, then for E > max{VC|_oo, V^-oo} one can 
set up a one-dimensional scattering problem in a completely standard manner — 
see, for example, standard references such as [3 [6], [71 [H [9], [101 [HI [I2l [13]. The 
scattering problem is completely characterized by the transmission and reflection 
amplitudes (denoted t and r), although the most important aspects of the physics 
can be extracted from the transmission and reflection probabilities (T = \t\^ and 
R= Irp). 



2.1 Ansatz 

The idea is to try to say things about exact solutions to the ODE 

i;"{x) + k'^{x)ij{x) = 0, (2.2) 

by comparing this ODE to some "simpler" one 

<(x) + e,ix) M^) = 0, (2.3) 

for which we are assumed to the know exact solutions ^^0(2;). In a manner similar to 
the analysis in references [H [2], we will start by introducing the ansatz 

ij{x) = a{x) i^oix) + b{x) ^poix). (2.4) 

This representation is of course extremely highly redundant, since one complex num- 
ber tlj{x) has been traded for two complex numbers a{x) and b{x). This redundancy 
allows us, without any loss of generality, to enforce one auxiliary constraint connect- 
ing a{x) and b{x). We find it particularly useful to enforce the auxiliary condition 

Subject to this auxiliary constraint on the derivatives of a{x) and b{x), the derivative 
of ip{x) takes on the especially simple form 

^ = a^l^'o + b ro'. (2.6) 

(This ansatz is largely inspired by the techniques of references [H [2], where JWKB 
estimates for the wave function were similarly used as a "basis" for formally writing 
down the exact solutions.) 
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2.2 Probability density and probability current 

For the probability density we have: 

p = ri' (2.7) 

= \a{x)tPo + Kx)ro\^ (2.8) 

= {\a\' + m\^|Jo\' + 2Re{ab*^|Jl} (2.9) 

= {|ar + |feri}po + 2Re{a6>o'}- (2-10) 

Furthermore, for the probabihty current: 



dx ^ 

= lm!^[a*ro + b*M [a< + bro']^ (2.12) 

= Im ||a|Vo>o + l^lVoV'o' + «&>o^o + a*6Vo>o'| (2.13) 

= {\a\'-mim{ro^'o} (2.14) 
= {\a\'-m^o. (2.15) 

Under the conditions we are interested in, (corresponding to a time-independent 
solution of the Schrodinger equation), we have p = 0, and so = 0. (And 

similarly po = 0, so d^^o = 0.) That is, ^ and are position-independent 
constants, which then puts a constraint on the amplitudes \a\ and \b\. Applying an 
appropriate boundary condition, which we can take to be a(— oo) = 1, b{—oo) = 0, 
we then see 

|ap-|6|2 = l. (2.16) 

This observation justifies interpreting a{x) and b{x) as "position-dependent Bogoli- 
ubov coefficients". Furthermore without any loss in generality we can choose the 
normalizations on ip and ipo so as to set the net fluxes to unity: ^ = = 1. 

2.3 Second derivatives of the wavefunction 

We shall now re-write the Schrodinger equation in terms of two coupled flrst-order 
differential equations for these position-dependent Bogoliubov coefficients a{x) and 
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b{x). To do this, evaluate d^^/^/da;^ making repeated use of the auxihary condition 

§ = ^(^^J + '-O (2.17) 

= a'^'o + b'^*' + a^'^ + b^*" (2.18) 

= a'i^'o-a'^ro' -akli^o-bklro (2-19) 
= ^ {^o>o - ^o^o'} - kl [a^o + b^o] (2.20) 



^0 



ki [a^Ao + bi\)l\ (2.21) 



2m' 

= —r-kl\a^^^HV\- (2.22) 
ro 

Where in the last line we have finally used our normalization choice = 1. This 
is one of the two relations we wish to establish. Now use the gauge condition to 
eliminate da/dx in favour of d6/dx to obtain a second relation for d'^tjj/dx'^. This 
now permits us to write d'^ip/dx'^ in either of the two equivalent forms 

dV 2za' , 2 



dx'^ ipQ 



-k',[atlJo + bro]; (2.23) 



2ib' 

■kl[aiJo + bro]- (2.24) 



^1 

2.4 SDE as a first-order system 

Now insert these formulae for the second derivative of the wavefunction into the 
Schrodinger equation written in the form 

d^^ 



+ k{xyip = 0, (2.25) 



dx'^ 

to deduce the pair of first-order ODEs: 

^ = +'^[k'-kl]{a\i;o\' + bro'h (2.26) 
£ = -'-[e-kl]{aijl + b\iJo\'}. (2.27) 

It is easy to verify that this first-order system is compatible with the auxiliary condi- 
tion (12. 5p . and that by iterating the system twice (subject to this auxiliary condition) 
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one recovers exactly the original Schrodinger equation. We can re-write this Ist-order 
system of ODEs in matrix form as 



dx 



ro' 1 



a 
b 



(2.28) 



(Matrix ODEs of this general form are often referred to as Shabhat-Zakharov or 
Zakharov-Shabat systems [1] . This matrix ODE can be used to write down a formal 
solution to the SDE in terms of "path-ordered exponentials" as in references [H |2] . 
We choose not to adopt this route here, instead opting for a more direct computation 
in terms of the magnitudes and phases of a and b.) 



2.5 Formal (partial) solution 

Define magnitudes and phases by 

a = \a\ e^"^"; b = \b\ e''^"; V'o = |^o| e^"^°. 

Calculate 



whence 



a' = \a\' e''^^ + i\a\ e'''''' (f)'^ = e'''''' {\a\' + i\a\ (p'J 



\a\' + t\a\ 01 = ^[k'- Q l^ol' {\a\ + |&| e-^(<^'^-^''+2^«)}. 



Similarly we also have 

\b\' + z\b\<P', 



a e 



-i(0(,-(/)a-2</)o)" 



Now take the real part of both these equations, whence 
, ,, 1 



[k-K] \b\ \^IJo\'sm{(l)a- 



Therefore 



That is 



= +2 - ^o] l«l l^ol' sin(0, - 06 + 20o). 
W\' = l[k^-kl] |^o|'sin(0,-0fe + 20o) V|a|2 - 1. 



[k -kl] I^Aol sin(0, -0f, + 20o), 



(2.29) 

(2.30) 

(2.31) 

(2.32) 

(2.33) 
(2.34) 

(2.35) 

(2.36) 



whence 



X 1 

{cosh~^ |a|}^^ = - / [k'^-k 

J X\ 



2] l^ol' sin(0, - 05 + 20o) dx. (2.37) 



Now apply the boundary conditions: At a; = — oo we have both a(— oo) = 1, and 
oo) = 0. Therefore 



cosh-^|a(x)| = ^ " 



[/c^ - /cq] IV^oT sin(0a - 0b + 20o) dx, 



;2.38) 



and so 



|a(x)| =cosh|^y" |^orsin(0,-0fe + 20o) dxj. (2.39) 



In particular 

cosh""*" |a(oo)| 
or equivalently 

lafcxo)! = cosh 



[/c^ - /c^] l^oP sin(0a - 06 + 20o) dx, 



[A;^ - Ajq] Iv^oT sin(0a - 06 + 20o) da; 



^2.40) 



(2.41) 



Of course this is only a formal solution since 0a (x) and 0;,(x) are, (at least at this 
stage), "unknown". But we shall argue that this formula still contains useful infor- 
mation. In particular, in view of the normalization conditions relating a and b, and 
the parity properties of cosh and sinh, we can also write 



|a(oo)| = cosh 



|fe((X))| = sinh 



- [k^- kl] IV^ol' sin(0a - 06 + 20o) dx 

^ ^-oo 

[fc^ - kl] l^oP sin (0a - 06 + 20o) dx 



■oo 
-oo 



oo 



(2.42) 
(2.43) 



2.6 First set of bounds 

To determine the first elementary set of bounds on a and h is now trivial. We just 
note that 

|sin(0a- 06 + 20o)| < 1. (2.44) 
Therefore ^ 

\a(oo)\ <coshi^ I \P - klWiJol^ dx\ ] (2.45) 



\b{oo)\ < sinh 1^ - kl\ IV^ordxj . (2.46) 

What does this now tell us about the Bogoliubov coefficients? 

2.7 Bogoliubov coefficients 

The slightly unusual thing, (compared to our earlier work in references [H [2], H]), is 
that now the "known" function tpo may also have its own Bogoliubov coefficients. 
Let us assume we have set our boundary conditions so that for the "known" situation 

ipoix ~ —oo) ~ exp{ik{—oo)x}, (2.47) 

and 

ipo{x ~ +oo) ~ ao exp{ik{+oo)x} + /3o exp{—ik{+oo)x} . (2.48) 
Then the way we have set things up, for the "full" problem we still have 

ipix ~ — oo) ~ exp{ik{—oo)x}, (2.49) 

whereas 

ip{x ^ +oo) ~ a{oo) ipoi^x) + b{oo) iPq{x) (2.50) 
~ [aoa{oo) + l3Qb{oo)] exp{ik{+oo)x} 

+ [/5o a(cxD) + 6(oo)] exp{—ik{+oo)x}. (2.51) 

That is, the overall Bogoliubov coefficients satisfy 

a = aoa(oo) +/5o 6(oo); (2.52) 

P = Poa{oo) + alb{oo). (2.53) 

These equations relate the Bogoliubov coefficients of the "full" problem {iplx), k{x)} 
to those of the simpler "known" problem {ipo{x), ko{x)}, plus the evolution of the 
a[x) and b{x) coefficients. Now observe that 

|a| < |ao| |a(oo)| + \f3o\ \b{oo)\. (2.54) 

But we can define 

|ao| = cosh9o; \Po\ = sinhBo; ^(oo)! = cosh 9; |&(oo)| = sinh 6; (2.55) 
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in terms of which 

|a| < cosh Oo cosh 6 + sinh 9o sinh 6 = cosh (Oq + 0) . (2.56) 
That is: Since we know 

/• + 00 

e < Gbound = 7^ \k^- kl\ \4Jo\^dx, (2.57) 



we can deduce 



(2.58) 



|a| < cosh |cosh ^ |ao| 2 J 1^^ ~ ^ol l^ol^da; j> ; 

<sinh|sinh-^|/3o| + ^y" lA;^ - |V^ordx| . (2.59) 
2.8 Second set of bounds 

A considerably trickier inequahty, now leading to a lower bound on the Bogoliubov 
coefficients, is obtained by considering what the phases would have to be to achieve 
as much destructive interference as possible. That implies 

|a| > |«o| |a(oo)| - \l3o\ Hoo)\, (2.60) 

whence 

|a| > coshlGo - e| . (2.61) 
Therefore, using 6 < 6bound, it follows that as long as 6bound < ©O; one can deduce 

|a| > cosh{eo-ebound}. (2.62) 

(If on the other hand Gbound > ©O; then one only obtains the trivial bound |a| > 1.) 
Another way of writing these bounds is as follows 

|a| > cosh Icosh""*^ |ao| ~ 2 J ~ ^ol iV^ol^dx j> ; (2.63) 

> sinh jsinh^i \(3o\ -\ f 1^^ " ^ol iV'oPdxj ; (2.64) 

with the tacit understanding that the bound remains valid only so long as argument 
of the hyperbolic function is positive. 



10 



2.9 Transmission probabilities 

As usual, the transmission probability (barrier penetration probability) is related to 
the Bogoliubov coefficient by 

2^ = (2-65) 



whence 



That is 



or even 



T > sech^ |cosh ^ |ao| ~^ 2 J 1^^ ~ ^ol IV^ordxj> 



(2.66) 



T > sech^ <j cosh~^(To"^/^) + 1 [ - kl\ iV^oPda; } , (2.67) 



2 



00 



T > sech^ |sech-i(ro/^) + \J^ 1^^ - kl\ \H^dxj . (2.68) 

Furthermore, as long as the argument of the sech is positive, we also have the upper 
bound 

T < sech^ |sech-^(To/') " ^ / " l'^' " '^ol IV^ordxj . (2.69) 

If one wishes to make the algebraic dependence on Tq clearer, by expanding the 
hyperbolic functions these formulae may be recast as 

To 



T > 



2 ' 



cosh || |P - Q l^ol'dx} + v^r^sinh {i - k^,\ IH'dxj 

(2>0) 

and (as long as the numerator is positive before squaring) 

To 



T < 



cosh {i |F - ki\ iV'olMx} - v/T^sinh {1 |F - ki\ l^olMx} 



2 • 



2.71 



3 Consistency check 

There is one special case in which we can easily compare with the previous results 
of references [H [2]. Take ko = k{±oo) to be independent of position, so that our 
comparison problem is a free particle. In that case 

^^^exp(p^. |V'or = ^; ^0 = 1; «o = l; Po = 0. (3.1) 
y/ko ko 
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Then the bounds derived above simphfy to 



|a|<cosh<j— / \k'^-kl\dx}, (3.2) 



<sinh<'-^-/ \k^-kl\dx}. (3.3) 



2k 



This is "Case I" of reference [T] and the "elementary bound" of reference [2], which 
demonstrates consistency whenever the formahsms overlap. (Note that it is not 
possible to obtain "Case II" of reference [1] or the "general bound" of reference P, |2] 
from the present analysis — this is not a problem, it is just an indication that this 
new bound really is a different bound that only partially overlaps with the previous 
results of references [H E], Hj. 

A second (elementary) check is to see what happens if we set ip{x) — ^ ipo{x)j 
effectively assuming that the full problem is analytically solvable. In that case T 
To, (and similarly both a — > ao and (3 — > Po), as indeed they should. 



4 Keeping the phases? 

We can extract a little more information by taking the imaginary parts of equations 
dMH) and fl^:^ to obtain: 



'l[k' - kl] l^or |l + 1^ cos(0, - 0, - 200 



(4.2) 



Subtracting 



(0a - 06)' = [k' - kl] |1 + cos(0, - 05 + 20o) | . (4.3) 

This is now a differential equation that only depends on the difference in the phases 
— the overall average phase (0a + 0fe)/2 has completely decoupled. (Likewise, in 
determining the transmission and reflection probabilities, this average phase also 
neatly decouples). To see how far we can push this observation, let us now define a 
"nett" phase 

A = 0a - 06 + 200. (4.4) 
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Furthermore, as per the previous subsections, we retain the definitions 

\a\ = cosh 6; \b\ = sinh 6. (4.5) 
Then equation f l2.39p becomes 

= U r - ^o] |V^orsin(A(x)) dx] . (4.6) 



while the "nett" phase satisfies 

We can even substitute for 0(x) and thus rewrite this as a single integro- differential 
equation for A(x): 

A{xy = {[e-kl] l^or +20'o} + . ^ cos[A(a:) 

sinh ( j_^[k'^ — k^] |-?/'oP sin[A(x)] dx 



(4.8) 

This equation is completely equivalent to the original Schrodinger equation we started 
from. Unfortunately further manipulations seem intractable, and it does not appear 
practicable to push these observations any further. 

5 Application: Small shift in the potential 

Let us now consider the situation 

V{x) = Vo{x) + e6V{x), (5.1) 

for e "sufficiently small" . 

5.1 First-order changes 

To be consistent with previous notation let us define 

e = kl + e j ^kl + e 6v. (5.2) 



Using equation (12.39^ we obtain the preliminary estimates 

|a(x)| = l + 0(e2), (5.3) 
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and similarly 

\b{x)\=0{e). (5.4) 

It is now useful to change variables by introducing some explicit phases so as to 
define 

a = d exp J[k'^ - k^] \iPq\ dx^ ; (5.5) 

b = bexp (^-'- J [k' - kl] l^o'l dxj . (5.6) 
Doing so modifies the system of differential equations (12.261 12.271) so that it becomes 

^ = +l[k' - kl] b exp {-I j - kl] dx^ ■ (5.7) 

' [k^ - kl] ~a ^l exp f +2 I \k^ - \^l\ dx^ . (5.8) 



dx 2 

The advantage of doing this is that in the current situation we can now estimate 



dx 
d6 

dx 2 



^ bv{x) i)l{x) exp ( +'i / e \iil\ dx \ + C(e^). (5.10) 



Integrating 



b[po) = J 6v{x) iplix) exp j e 5v \%l)l\ dx^ dx + C(e^). (5.11) 

This is not the standard Born approximation, though it can be viewed as an instance 
of the so-called "distorted Born wave approximation" . In terms of the absolute values 
we definitely have 

|6(oo)| = |6(oo)| < - / \5v{x)\\^l{x)\dx + 0{e'^). (5.12) 



oo 



5.2 Particle production 

When it comes to considering particle production we note that 

/? = /3o a(oo) + al b{oo) = f3o + a* b{oo) + 0{e^), (5.13) 
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so the change in the number of particles produced is 

5\(5^\ = Re{2a*/3oKoo)} + C(e^). (5.14) 

In particular 

/ + 00 

\5v{x)\\^l{x)\dx + 0{e^), (5.15) 
-co 

which we can also write as 

r+oo 

\5N\<e v/iVo(A^o + l) / \5v{x)\\^Pl{x)\ + 0{e^). (5.16) 



Note that one will only get an order e change in the particle production if the "known" 
problem {ipQ, k^} already results in nonzero particle production. 

5.3 Transmission probability 

To see how a small shift in the potential affects the transmission probability we note 
But then 

T = 2, (5.18) 

|aoP |l + {/?oKoo)/«o} + 0(e2)r 



implying 



So the change in the transmission probability is 



T = To <{ 1 - 2Re <! K + o{e') \ . (5.19) 



5T = -To i 2Re <{ ^2-^ K + O(e^) } . (5.20) 



"0 

Taking absolute values one obtains 

/ + 00 
\6vix)\ \^l{x)\dx + 0{e'). (5.21) 
-oo 

Note that one will only get an order e change in the transmission probability if 
the "known" problem {ipo, ko} already results in nonzero transmission (and nonzero 
reflection). 
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6 Discussion 



What are the advantages of the particular bounds derived in this article? 

• They are very simple to derive — the algebra is a lot less complicated than 
some of the other approaches that have been developed P, El [3l H] . (And a lot 
less complicated than some of the blind alleys we have explored.) 

• Under suitable circumstances the procedure of this article yields both upper 
and lower bounds. Obtaining both upper and lower bounds is in general very 
difficult to do — see in particular the attempts in [2] . 

• All of the other bounds we have developed [H E], [3], Hj needed some condition on 
the phase of the wave-function, (some condition similar to 7^ 0), which had 
the ultimate effect of making it difficult to make statements about tunnelling 
"under the barrier". There is no such requirement in the present analysis. 
(The closest analogue is that we need 7^ 0, which we normalize without 
loss of generality to ^0 = !•) In particular this means that there should be no 
particular difficulty in applying the bound in the classically forbidden region 
— the "art" will lie in finding a suitable form for t\}Q which is simple enough to 
carry out exact computations while still providing useful information. 

In closing, we reiterate the fact that generic one-dimensional scattering problems, 
which have been extensively studied for close to a century, nevertheless still lead to 
interesting features and novel results. 
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